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A  model  has  been  developed  that  can  be  used  to  predict  the  polarization 
state  of  scattered  light  from  surface  particles  satisfying  the  specific 
criterion  that  the  polarization  state  of  the  light  scattered  by  the  isolated 
particle  is  symmetric  and  resembles  that  of  a  Rayleigh  sphere.  Although 
this  condition  appears  constraining,  it  may  be  satisfied  by  many  common 
contaminants  of  interest,  such  as  dust  and  aggregates  of  Rayleigh-sized 
particles.  When  these  contaminants  are  placed  on  a  substrate  and  the 
system  is  illuminated  at  near-grazing  incidence,  the  polarization  state  of 
the  scattered  light  is  shown  to  be  nearly  independent  of  the  actual  shape 
and  refractive  index  of  the  contaminant.  It  is,  however,  sensitive  to  the 
dimensions  of  the  particle  and  substrate  refractive  index.  Thus,  the  model 
may  be  useful  in  determining  the  size  of  surface  contaminants. 
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1.  Introduction 


\ 


Many  contaminants  of  interest  are  irregular  particles,  like  dust  and  soot, 
which  might  be  modeled  as  aggregates.  Smoke  particles,  for  instance, 
coagulate  into  sparse  random  clusters  of  Rayleigh-sized  spherules. 
Although  such  clusters  are  large  compared  to  the  wavelength,  their 
polarization  Mueller  matrix  elements  (S12,  S33,  and  534)  resemble  those  of 
a  Rayleigh  particle  [1].  This  similarity  is  also  seen  for  complex  structures 
like  Arizona  road  dust,  a  sample  of  which  is  shown  in  figure  1.  Mueller 
matrix  elements  for  Arizona  road  dust  have  been  measured  with  the 
polarization  modulation  techniques  developed  by  Hunt  and  Huffman  [2] 
and  are  compared  to  the  elements  of  a  Rayleigh  sphere  in  figure  2.  Even 
though  the  dust  particles  are  many  times  larger  than  the  wavelength,  the 
polarization  matrix  elements  have  characteristics  similar  to  those  of  the 
Rayleigh  sphere:  the  characteristic  dip  in  matrix  element  S12,  the 
characteristic  S-shape  of  element  533,  and  element  £34  ~  0.  The  evidence 
of  the  particle's  size  is  contained  within  the  total  intensity  element  Sn. 
Rather  than  a  symmetric  intensity  distribution  with  only  a  50-percent 
dropoff  at  6  =  90°  (as  would  be  expected  for  a  Rayleigh  sphere),  we  see  a 
three-order-of-magmtude  drop  in  intensity  from  the  forward  scatter, 
followed  by  an  order-of-magnitude  rise  in  intensity  in  the  backscatter. 
Even  though  the  total  scattered  intensity  does  not  resemble  that  of  a 
Rayleigh  particle,  the  polarization  elements  strongly  mimic  Rayleigh 
scattering.  I  use  this  observation  in  developing  a  model  of  the  polarization 
state  of  the  light  scattered  from  such  particles  placed  on  a  substrate.  Since 
irregular  particles  like  dust  and  aggregates  are  common  byproducts  in 
many  manufacturing  processes,  characterizing  such  particles  is  one  of  the 
motivating  factors  in  developing  such  scattering  models  [3]. 
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Figure  1.  Arizona  road 
dust. 


2.  Model 


\ 


One  method  of  solving  for  the  scattered  fields  from  a  particle  near  a  plane 
interface  is  to  use  image  theory.  Figure  3  shows  the  geometry  of  the 
scattering  system.  A  particle  is  placed  a  distance  d  beneath  a  plane 
interface  separating  a  substrate  of  refractive  index  nsub  (above)  from  air 
(below).  The  incident,  scattered,  and  interaction  fields  are  each  expressed 
as  multipole  expansions  of  vector  spherical  harmonics,  and  the  boundary 
conditions  at  the  surfaces  of  the  particle  and  of  the  substrate  are  then 
satisfied  simultaneously  [4-8]. 

Consider  expansions  about  the  particle  coordinate  system  (j  =  1)  and  the 
image  coordinate  system  (j  =  2).  These  fields  expanded  about  the  particle 
coordinate  system  are 


Ew  =  V  V  M(1)  ,  +  a<2>  N(1)  , 

(i) 

n=0  m——n 

oo  n 

E*»  =  E  E 

(2) 

71=0  TTl— — 71 

oo  n 

Eint  =  V  V  c^1)  M(1)  ,  +  CW  N(1)  , 

(3) 

n=0  m=—n 


where  and  j  are  vector  spherical  harmonics  expanded  by 

Bessel  functions  of  the  ith  kind  expanded  about  the  j  =  1  coordinate 
system.  The  scattering  solution  can  be  found  directly  from  the  scattering 
response,  or  T-matrix,  of  the  isolated  particle  to  an  incident  wave  of  the 
form 

Einc  =  4^M^  +  429)N$  (4) 

and  is  given  by 


E 


sea  — 


E 

n',771' 


(5) 


J 


The  solution  to  the  scatter  from  an  arbitrarily  shaped  particle  placed  near 
a  plane  interface  has  recently  been  formulated  [7,8].  These  results  can 
easily  be  modified  to  accommodate  substrates  that  are  not  perfectly 
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Figure  3.  Geometry  of 
scattering  system,  showing 
particle  placed  below 
substrate  centered  on 
coordinate  system  j  =  1.  Its 
image  is  centered  above 
substrate  on  coordinate 
system  j  =  2. 


conducting,  with  the  following  approximation:  say  the  interaction  field  is 
the  image  of  the  scattered  field  from  the  particle  that  strikes  the  substrate 
at  normal  incidence  before  again  interacting  with  the  particle  [4,6].  The 
interaction  coefficients  can  then  be  expressed  as 


cW  =  Rte  (0)  £  (-1  )”'+m  -  »SXn>)] .  («) 

n' 

=  Rte  (0)  £  (-1  )”,+m  ■ ”>] ,  (7) 

n' 

where  D^1'^  and  are  translation  coefficients  used  to  translate 

vector  spherical  harmonics  in  the  image  coordinate  system  to  vector 
spherical  harmonics  in  the  particle  coordinate  system: 
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and  Rt*  {0)  are  the  standard  Fresnel  reflection  coefficients  for  a  plane 
wave  reflecting  off  the  substrate.  The  scattering  coefficients  are  then 


(9) 
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Adding  the  electric-field  components  emanating  from  the  particle  and  its 
image  yields  the  scattering  amplitude  matrix: 

E#00  \  _  exp[*Ax] 
j 5«ca  I  —  ikr 

where  the  scattered  amplitudes  are  the  superposition  of  the  scattered 
fields  from  the  object  and  the  image: 
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The  incident  field  coefficients  bnmE  and  bnmM  in  these  expressions  are 
for  the  superposition  of  the  real  and  image  TE  and  TM  plane  waves  [4]. 
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3.  Some  Approximations 


The  equations  in  the  previous  section  may  be  used  to  calculate  the 
scattered  fields  from  a  particle  near  a  highly  reflecting  substrate.  Their 
primary  drawback  is  that  numerically  finding  the  scattering  coefficients 
(eq  (10))  requires  a  matrix  inversion.  Although  such  numerical  calculations 
are  achieved  more  cost-efficiently  than  are  numerical  techniques  like 
discrete  dipole  approximation  (DDA),  they  are  still  not  conducive  to  the 
rapid  calculations  necessary  for  modeling  surfaces  produced  in  the 
manufacturing  process.  More  importantly,  when  such  operations  are 
performed,  much  of  the  physical  insight  into  a  solution  is  lost. 

In  certain  instances,  the  matrix  is  diagonal,  and  the  solution  is  greatly 
simplified.  If,  for  instance,  the  interaction  field  is  small  compared  to  the 
incident  field,  the  off-diagonal  elements  of  the  matrix  form  of  equation 
(10)  are  small  compared  with  the  diagonal  elements,  and  bnm  ~ 

Such  conditions  obviously  exist  for  Rayleigh  and 
Rayleigh-Gans  particles  and  particles  placed  a  large  distance  from  the 
substrate.  The  scattered  field  from  a  particle  on  a  substrate  can  then  be 
found  directly  as  the  superposition  of  the  scattered  fields  from  the  particle 
and  the  image  particle,  without  any  consideration  of  the  interaction.  In  the 
plane  defined  by  the  surface  normal  and  the  incident  wavevector,  the 
scattered  field  from  a  particle  resting  on  a  substrate  can  be  expressed  as 
the  superposition  of  the  fields  from  the  particle  and  image  particle;  i.e.,  in 
the  plane  of  incidence,  equations  (12)  and  (13)  may  be  written  as 


51  (0  -  0O)  =  S[  ( 0O ,  9  -0O)  +  S°  (tt  -0O,9  +  9O-  n r)  Rte  (#o)  exp  [iJ  (0O,  d )]  (16) 

+  S°  ( 6a ,  7 t  -0  -0o)  Rte  (tt  -  0)  exp  (7 r  -  0,  d)] 

+  S°  (ir  -  60, 0o  -  0)  RTe  (#<>)  Rte  (it  -  0)  exp  {i  [<5  ( 0O ,  d)  +  6  (n  -  0,  d)]}  , 

52  ( 0  -  90)  =  S%  (00, 0-9o)  +  S°2  (7 r  -  9a,  9  +  9a  -  7 r)  RTM  (#0)  exp  [iJ  (0O,  d)]  (17) 

+  S2  (9o,  7 x  —  0  —  0O)  Rtm  (tt  —  9)  exp  [i5  (-n  —  9 ,  d)] 

+  S2(n-  00, 0O  -  9)  Rtm  (0»)  Rtm  (*  ~  0)  exp  {i  [<5  ( 0O ,  d)  +  5  (n  -  0,  d)]}  , 

where  in  the  far  field  9  ~  9\  ~  02,  6  (9,  d)  =  2kd  cos  9)  also.  Si  (00, 9  —  0O) 
and  £2  (6*0, 0  -  0O)  are  the  scattering  amplitude  matrix  elements  for  the 
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isolated  particle  illuminated  by  a  plane  wave  whose  wavevector  is 
oriented  at  angle  0O  with  respect  to  the  z-axis,  and  9  —  0O  is  the  scattering 
angle  measured  from  the  specular  direction.  Equations  (16)  and  (17)  are  a 
direct  result  of  setting  the  interaction  fields  equal  to  zero  in  the  equations 
describing  the  scatter  from  particles  on  or  near  substrates  [4-12]. 

It  is  likely  that  the  interaction  fields  will  be  negligible  for  some  highly 
irregular  particles  resting  on  a  substrate:  e.g.,  particles  consisting  of 
aggregates  of  Rayleigh  spheres.  The  polarization  state  of  the  scattered 
light  from  these  isolated  aggregates  resembles  that  of  a  Rayleigh  sphere; 
i.e.,  the  particle  acts  as  a  group  of  individual,  noninteracting  Rayleigh 
subparticles,  for  which  there  is  interference  between  the  waves  emanating 
from  each  subparticle  but  very  little  interaction  between  the  subparticles 
[1].  Since  there  is  negligible  interaction  between  subparticles  of  the 
system,  it  can  reasonably  be  assumed  that  there  would  also  be  negligible 
interaction  between  the  particle  and  its  image  when  it  is  placed  on  a 
substrate.  The  scattering  amplitude  matrix  elements  for  an  isolated 
Rayleigh  scattering  system  are  related  by 

S%  (0O,  $-0o)=  cos  (0  -  0O)  5?  (00, 0  -  0O)  •  (18) 

This  equation  is  also  valid  for  a  conglomeration  of  noninteracting 
Rayleigh  spheres  and  is  approximately  valid  for  some  highly  irregular 
particles  (e.g.,  the  Arizona  road  dust  of  fig.  1  and  2).  By  using  the 
Rayleigh-scattering  approximation  of  equation  (18),  one  can  express 
S%  (00, 0  —  0o)  in  terms  of  Sf  ( 00, 0  —  0O)  in  equations  (16)  and  (17). 

One  final  condition  in  developing  this  scattering  model  is  that  at  some 
particle-coordinate-system  location,  the  scattering  amplitude  matrix 
elements  for  the  isolated  particle  are  symmetric: 

s°i(do,o-o0)  =  s°1(e0,e0-e).  (is>) 


Such  a  condition  is  valid  for  symmetric  particles  such  as  spheres, 
ellipsoids,  and  dipoles  when  the  coordinate  system  is  at  the  center  of  the 
particle.  In  addition,  we  know  that  the  scattering  phase  function  (which  is 
the  sum  of  the  magnitudes  squared  of  the  scattering  amplitude  matrix 
elements)  of  many  particles  is  also  approximately  symmetric.  The 
scattering  amplitude  matrix  elements,  however,  contain  the  phase 
information  of  the  scattered  fields  and  are  dependent  on  the  placement  of 
the  particle  coordinate  system.  For  example,  if  we  move  the  coordinate 
system  close  to  one  edge  of  a  spherical  scatterer,  the  complex  amplitudes 
of  the  scattered  spherical  waves  expressed  in  that  coordinate  system  are 
no  longer  symmetric,  even  though  the  intensity  distribution  in  the  far  field 
is  quite  symmetric. 
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For  highly  irregular  systems,  equation  (19)  may  hold  some  validity  at  a 
coordinate  system  placed  at  some  central  location  of  the  particle.  The 
accuracy  of  using  this  approximation  is  further  increased  if  we  average  the 
scattering  amplitudes  over  some  finite  area  to  average  out  the 
high-frequency  interference  structure  contained  in  the  scattering 
amplitude  coefficients.  This  is  physically  equivalent  to  having  collection 
optics  or  a  finite-size  detector  collect  the  scattered  field.  Although  not  all 
irregular  particles  would  be  expected  to  satisfy  this  condition,  the 
comparison  with  experimental  results  given  below  shows  that  this 
condition  is  at  least  partially  satisfied  for  some  substrate  contaminants. 

We  can  now  derive  the  polarization  state  of  scattered  light  for 
contaminants  fulfilling  the  symmetric-scatterer  (eq  (19))  and  the 
noninteracting-Rayleigh-scatterer  (eq  (18))  conditions.  If  the  substrate  is 
illuminated  at  near-grazing  incidence  ( 0o  ~  ±90°),  the  scattering 
amplitudes  of  the  isolated  particle  (S°(0O, 6  -  0o)  and  S%(0o,  0  -  0o))  may 
be  factored  out  of  equations  (16)  and  (17),  and  completely  cancel  from  the 
normalized  polarization  matrix  elements.  The  normalized  polarization 
matrix  elements  now  depend  only  on  the  separation  distance  d  of  the 
particle  from  the  substrate,  the  substrate  complex  refractive  index  nsub, 
and  the  incident  angle  0o  (as  long  as  6a  ~  ±90°). 

It  is  important  to  note  that  the  separation  parameter  d  is  the  distance  from 
the  surface  of  the  substrate  to  some  particle  coordinate  system  location  at 
which  the  symmetric-scatterer  approximation  holds.  The  physical  extent 
of  a  symmetric  scatterer  is  equal  to  2 d,  since  this  location  must  be  at  the 
center  of  such  a  particle.  If  the  particle  is  not  symmetric,  we  would  still 
expect  it  to  be  located  near  some  central  location  of  the  particle,  so  that  the 
particle  would  have  a  physical  dimension  of  approximately  2d.  With  a 
particle  of  this  size,  we  can  now  apply  the  model  to  determine  the  size  of 
contaminants  located  on  substrates  by  varying  the  separation  distance 
(particle  size)  in  the  model  until  the  polarization  state  of  the  scattered  light 
matches  the  experimentally  measured  state.  Finally,  note  that  illuminating 
the  system  at  grazing  incidence  also  reduces  the  interaction  term,  since  the 
large  forward-scattering  lobe  is  directed  along  the  surface  rather  than  back 
toward  the  particle. 

To  model  the  scattered  light  of  contaminants  on  a  substrate,  let  us 
compare  model  results  with  the  normalized  polarization  Mueller  matrix 
elements  of  select  perfect  and  decayed  optical  substrates  measured  by 
Iafelice  and  Bickel  [13].  These  systems  are  Al-coated  sapphire  substrates 
with  rms  surface  roughness  a  =  1.6  ±  0.3  nm.  Flowever,  the  authors  note 
that  surface  defects  consisting  of  residual  particles  ranging  in  size  from 
~0.2  to  1.0  /xm  sparsely  littered  both  substrates.  Although  no  difference  in 


size  distribution  was  apparent,  the  decayed  mirror  surface  had 
approximately  five  times  as  many  defects  as  the  perfect  mirror  surface. 
During  the  measurements,  the  laser  light  was  focussed  at  several  different 
locations  on  the  substrate,  and  sites  with  large  scattering  were  avoided; 
thus,  it  is  likely  that  only  one  of  the  smaller  scatterers  was  illuminated 
when  measurements  were  taken.  Figure  4  shows  the  Mueller  matrix 
elements  measured  [13]  at  near-grazing  incidence  ( 0o  =  79°)  with  A  = 
0.4416  /im.  Also  shown  are  the  calculated  polarization  matrix  elements  for 
a  symmetric  scatterer  at  a  separation  distance  d  =  A/4.7  from  the  (nsub  = 
0.5  +  3.5?)  substrate.  Since  the  substrate  refractive  index  is  highly 
dependent  on  the  thickness  of  the  aluminum  oxide  layer  that  forms  on  the 
surface  and  acts  as  an  anti-reflection  coating,  this  value  is  reasonable. 
Values  for  pure  A1  are  on  the  order  of  nsub  =  0.5  +  5.0?.  Videen  and  Ngo 
[14]  compared  theoretical  and  experimental  Mueller  matrix  elements  of  a 
cylinder  resting  on  a  substrate  prepared  by  the  same  method  as  Iafelice 
and  Bickel  and  found  best  fits  for  nsub  =  0.5  +  3.8?.  The  separation  distance 
(d  =  A/ 4.7)  corresponds  to  particles  having  dimensions  of  ~0.2  //of  the 
estimated  size  range.  Overall,  the  fit  of  the  model  results  to  measurement 
is  good.  Model  calculations  of  matrix  element  S12  bear  a  striking 
resemblance  to  those  of  the  decayed  mirror  surface  (increasing  the 
separation  distance  slightly  also  produces  the  hump  at  0  -  0o  ~  85°  for  the 
perfect  mirror  surface).  The  general  shapes  of  measurement  and  model 
results  for  matrix  elements  S33  and  S34  are  also  in  good  agreement. 
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Figure  4.  Normalized 
polarization  Mueller 
matrix  elements  calculated 
for  a  symmetric,  irregular 
scatterer  at  d  —  A/4.7  from 
a  substrate  (nsub  =  0.5  - 
3.5i)  illuminated  at  0o  = 
79°  (dots).  Superimposed 
are  experimental  matrix 
elements  measured  by 
Iafelice  and  Bickel  [13]  of  a 
perfect  mirror  surface 
(solid)  and  a  decayed 
mirror  surface  (dashes), 

A  =  0.4416  /zm. 
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4.  Conclusion 


A  model  has  been  developed  that  can  be  used  to  predict  the  polarization 
state  of  light  scattered  from  irregular  particles  on  a  plane  interface 
illuminated  at  near-grazing  incidence.  The  model  assumes  that  the 
interaction  of  the  particle  with  its  image  is  small  compared  with  the 
illuminating  field  and  that  the  scattered  field  from  the  isolated  particle  is 
symmetric  in  some  particulate  coordinate  system.  The  polarization  state  of 
the  scattered  light  is  sensitive  to  the  particle  size  (distance  from  the 
substrate  to  the  particulate  coordinate  system)  and  the  substrate  refractive 
index,  but  is  insensitive  to  the  actual  shape  of  the  contaminant.  Such  a 
model  is  useful  for  sizing  contaminants. 
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